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Abstract 

£ 

We investigate a non solvable two-dimensional ferromagnetic Ising 
model with nearest neighbor plus weak finite range interactions of 
j>. strength A. We rigorously establish one of the predictions of Confor- 

mal Field Theory (CFT), namely the fact that at the critical temper- 
ature the finite size corrections to the free energy are universal, in the 

^r) sense that they are exactly independent of the interaction. The corre- 

sponding central charge, defined in terms of the coefficient of the first 

t-h subleading term to the free energy, as proposed by Affleck and Blote- 

Cardy-Nightingale, is constant and equal to 1/2 for all < A < Ao and 
Ao a small but finite convergence radius. This is one of the very few 

. ^ cases where the predictions of CFT can be rigorously verified starting 

^ from a microscopic non solvable statistical model. The proof uses a 

combination of rigorous renormalization group methods with a novel 
partition function inequality, valid for ferromagnetic interactions. 
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1 Introduction and main results 

The applications of Conformal Field Theory (CFT) to statistical mechanics 
are based on the assumption that a statistical model at the critical point 
admits a non-trivial, conformally invariant, scaling limit, as suggested by the 
renormalization theory of critical phenomena. The two-dimensional (2D) lo- 
cal scale invariance strongly constraints the structure of the critical theory, 
as understood by Belavin, Polyakov and Zamolodchikov [2] . They recognized 
that the theory is characterized by a dimensionless constant c, the central 
charge, which is associated with an anomaly term in the commutation rela- 
tions of the stress energy tensor. The central charge can be also defined in 
terms of the finite size corrections to the free energy at criticality [HE]- In 
some cases, the critical theory is fully characterized by the value of c, which 
takes the form c = 1 — G/m(m + 1), m = 2, 3, 4, . . . Once c is fixed to one of 
these special values, the critical exponents are all explicitly known in terms 
of the Kac formula [TTJ. 

In practice, the identification of the critical theory associated with a given 
microscopic lattice model is done by inspection, by trying to match the known 
informations about the lattice model's exponents with the Kac formula. Once 
a correspondence is established or guessed, a large number of non trivial 
predictions on the model's correlation functions at criticality can be inferred, 
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which in general cannot be analytically derived by other means. It is therefore 
important to check these predictions in specific models, which could serve as 
benchmarks for this scheme. Unfortunately, there are just a few cases, based 
on exactly solvable lattice models, in which so far this correspondence could 
be rigorously established. A remarkable example is the nearest neighbor Ising 
model at the critical point, whose understanding in the scaling limit improved 
substantially in the last few months [7J El [9] . More in general, it is very hard 
to rigorously compute critical exponents, correlation functions or finite size 
corrections to thermodynamic functions at the critical point. In recent times 
new methods for the analysis of non-integrable 2D spin systems, based on 
the Renormalization Group, have been developed, starting from [21] and [2T] , 
where the authors computed the critical exponent of the so-called energy field 
operator in a class of perturbed Ising models (such as the one considered in 
this paper) and in a class of two stacked interacting Ising models (including 
the 8 vertex and the Ashkin- Teller models), respectively. By these methods 
one can try to verify some of the CFT predictions in the context of non 
solvable lattice models. 

In particular in this paper we consider an Ising model with a generic 
ferromagnetic short range interaction, with Hamiltonian 

H = - J ^ ^2 ^x+ej ~ * °x w ( x ~ y)°y ' ( L1 ) 

xGA £ , L j=l,2 {x,y} 

where J and A are positive constants, A^l C Z 2 is a finite rectangular box of 
sides L and £ with periodic boundary conditions, cr x G {±1}, and e,- are the 



two unit coordinate vectors on Z 2 . The sum in the second term of Eq.(l.l) is 
over all unordered pairs of sites in A^; the interaction potential v(x — y) is 
rotation invariant, positive and has finite range, namely: u(x— y) = 0, V|x| > 
Ro '■— Mo, for a suitable positive integer M . With no loss of generality, we 
can assume that v(0) = v(ej) = 0. The case A = corresponds to the 
nearest-neighbor Ising model which is exactly solvable [221 EE 1201 [23 122]; in 
the case A ^ no solution is known but in [21] the exponent of the energy- 
energy correlation was computed and shown to be universal (i.e., equal to 
2 as in the A = case), in contrast with the critical temperature or the 
amplitude of the correlations which are model dependent. 

A number of key informations on the system are encoded in the partition 
function: 

Zp(A e>L ) = e ~ m<L) ■ (1-2) 
At all temperatures, the thermodynamic limit for the pressure is well-defined 
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and independent of the speed at which £ and L are sent to infinity: 



l,L— >oo lL 

While the mere existence of the limit can be proved very generally, based 
on convexity and subadditivity arguments [25], the explicit form of pp can 
be computed either by the Onsager's solution (for A = 0) or by cluster 
expansion and Renormalization Group analysis (for A sufficiently small, see 
below). It turns out that pp depends explicitly on A, at all temperatures 
p. Moreover, at all temperatures but the critical one, the limit is reached 
exponentially fast. At f3 = /3 C (A) the limit is reached polynomially, and we 
shall define := Ppj\)- Remarkably, according to the ideas and methods 
of CFT, the finite size corrections to /oo are expected to be universal, in 
particular independent of A and v. More precisely, in the presence of periodic 
boundary conditions, Ref-JUE] predicted the validity of the following formula, 
asymptotically for large I: 

tZ II = f » + + V • (L4) 

where c is the central charge of the critical theory. For the nearest neighbor 
Ising model c = 1/2 and supposedly the same should be true for perturbed 



Ising models of the form (1.1). Our main results is a rigorous confirmation 
of this expectation. 



Theorem 1.1. Given the model (1.1) with A positive and small enough, there 



exists a critical temperature /3 C (A) such that Eq.(1.4) holds with 



c = l/2 (1.5) 

Recall that the critical exponent of the energy field operator is known 
to be universal and equal to 2 (2U E], as expected from the use of the 
Kac formula at c = 1/2 (see e.g. [TT]): therefore, our result says not only 
that the finite size corrections to the free energy are universal at criticality 
and that the corresponding central charge is constant, but also that such 
central charge has the right value, i.e. the one matching with the one guessed 
from the critical exponents and the Kac formula. This is a very non trivial 
connection, predicted by CFT, between critical exponents and finite size 
corrections to the thermodynamic functions, and our results is the first proof 
of the correctness of this prediction in a non-integrable statistical model. It 
would be very interesting to extend the connection to the critical exponent 
of the spin field operator, but this seems a much harder problem than the 
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one solved here; on the basis of the above correspondence, such exponent is 
expected to be 1/4, but its rigorous computation at A 7^ is still beyond 
reach of the current techniques. Another very interesting extension would be 
to prove an analogous universality result for statistical models with c = 1, 
such as the 8 vertex or the Ashkin- Teller model. In these cases several critical 
exponents have been computed and proved to be A-dependent [211 HHJ IH E] , 
as expected from CFT; at the same time, the finite size correction to the free 



energy are expected to be independent of A and given by Eq.(1.4) with c = 1, 
but this fact is unproved so far. We hope to come back to this question in a 
future publication. 



In the case of the nearest-neighbor Ising model, A = 0, the result Eq.(1.5) 
was proved by [10] and follows from the exact solution. If A 7^ 0, the proof 
of our main theorem uses a Renormalization Group (RG) analysis first in- 
troduced in [21] for the computation of the two-point energy correlation 
function, and recently extended in [H] to the analysis of the n-point energy 
correlations. In addition to the ideas of [2H H3], we use here some novel 
partition function inequalities, which we can only prove for ferromagnetic in- 
teractions, A > 0. The assumption of positivity of the interaction is expected 



to be technical and we believe that the analogue of Theorem 1.1 should hold 
also for A small and negative. 

More in detail, our proof proceeds as follows. We start by writing the 
partition function of the interacting Ising model as a sum of four Grass- 
mann integrals with different boundary conditions, in a way that naturally 
extends the analogous representation for the nearest neighbor model. Us- 
ing a partition function inequality, we reduce ourselves to the study of just 
one out of these four Grassmann partition functions, namely the one with 
antiperiodic boundary conditions; we prove that the other three terms are 
subleading and do not contribute to the central charge. On the other side, the 
Grassmann partition function with antiperiodic boundary conditions can be 
rewritten as the product of two terms: the first is equal to the non-interacting 
Grassmann partition function with renormalized parameters (this is the con- 
tribution from the "infrared fixed point"), while the rest includes all the 
corrections coming from finite infrared scales and from the irrelevant terms. 
Now, remarkably: 

1. the dependence of the first factor upon the renormalized parameters can 
be scaled out by a simple change of variables, after which the factor 
takes the form z 2xVolume x (free partition function), with Z the "wave 
function renormalization", which is volume-independent; therefore, the 
presence of the renormalized parameters changes the bulk pressure, but 
not the finite size corrections; 
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2. the rest can be studied by a multiscale analysis, which requires the 
introduction of two running coupling constants, playing the role of wave 
function and critical temperature counterterms; these running coupling 
constant go to zero exponentially fast in the infrared limit, thanks to a 
dimensional improvement in the dimensional bounds following from the 
fact that the theory is super-renormalizable, in the sense that all the 
field operators with more than two fields are irrelevant in the RG sense; 
correspondingly, the finite size corrections to the pressure coming from 
these term can be shown to go to zero faster than i~ 2 as the infrared 
cutoff I is removed, namely like £~ 2 ~ e , for < 9 < 1. 



The strategy resembles closely the one used in [TBI [T7] to prove the univer- 
sality of the optical conductivity in interacting graphene. An extra difficulty 
that we have to face in our case is the definition of the localization and renor- 
malization procedure at finite volume, which uses and extends the strategy 
proposed in j3]. 

The rest of the paper is organized as follows. In Section [2] we review the 
proof of Theorem 1.1 in the non-interacting case A = 0. In Section [3] we prove 



the main theorem in the A ^ case: we first review the Grassmann repre- 



sentation of the interacting partition function (Section 3.1); then we state 
the partition function inequality that we use to effectively eliminate three 



out of the four Grassmann partition function (Section 3.2); next we describe 
the RG computation of the antiperiodic Grassmann partition function (Sec- 



tion 3.3), and we use it to compute the bulk and subleading contributions 
to the pressure (Section 3.4). In Section [i] we prove the partition function 
inequality stated in Section |3.2| In the Appendix, we study a subleading 



correction to the pressure coming from the ratio of the Grassmann partition 



functions, both in the non-interacting (Appendix A.l) and in the interacting 



case (Appendix 3.4). 



2 The nearest-neighbor Ising model 

In this section we review the proof of our main theorem in the case of the 
nearest neighbor (n.n.) Ising model, A = 0. The proof can be found in [10] , 
but is reproduced here for the reader's convenience. From now on we shall 
drop the dependence on (3 in the symbol used for the partition function, 
in order to avoid a too cumbersome notation. We shall be as explicit as 
possible in distinguishing the formulas where the value of (3 is generic from 
those where (3 is fixed to be the critical one. We shall also assume that i and 
L are even, in order to simplify the signs appearing in some formulas. 
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The starting point is the representation of the n.n. Ising model's par- 
tition function Z°({ J b }\ A^i) with periodic boundary conditions and bond- 
dependent link variables in terms of a sum over multipolygons (see e.g. [T2]): 

Z\{J b };A^ L ) = 2 eL [ J] cosh(M)] yi niI tanh ^ Jb )' 

where B^l is the set of n.n. bonds in A^x and T is a collection of disjoint 
polygons in A^l, each of which is a closed connected collection of bonds; 
here closed means that at each of the vertices covered by the polygon there 
is an even number (either 2 or 4) of incident bonds. Note that some of the 
polygons in T can wind up the torus A^l, due to the periodic boundary 
conditions. Note also that the case T = is included in the sum, in which 
case the corresponding contribution is equal to 1. 

A convenient (for computational purposes) way of re-expressing the par- 



tition sum Eq.(2.1) is by writing it in terms of Grassmann integrals, see e.g. 



EaUH. Define 

Z° a ({Jb};A e , L ) = 2* L [ J] cosh(M)] [v<S>e s ^, 

S({J b };<5>) = J2 [tanh(/3J ( x,x+ei))-^x-£^x+ei + tanll(/3 J( X; x+e 2 ) ) ^x ^x+e 2 



xeA 

+HA + V X V X + V X H X + V X H X + H X V X + V X H X 



(2.2) 



Here H x , H x , V x , V x are independent Grassmann variables, four for each lat- 
tice site, <3> = {H x , H x , V x , ^4;} x eA denotes the collection of all of these 
Grassmann symbols and P<3? is a shorthand for f| x dH x dH x dV x dV x . The 
label a. = (0:1,02), with 01,02 £ refers to the boundary conditions, 

which are periodic or antiperiodic in the horizontal (resp. vertical) direction, 
depending on whether 01 (resp. o 2 ) is equal to + or — . The connection 
between the Grassmann integral and the mutlipolygon representation can be 
made apparent by expanding the exponential inside the Grassmann integral 
and by integrating term by term. The result is 

Z° a ({J b }; A e<L ) = 2 eL [ H cosM/U)] ■ (2.3) 

b£B e ,L 

• j2 (-ai)^(-^(-i) h ^nn tanh ^), 

rcA fjL 7er be-y 

where h(T) and v(T) are the number of windings of T on the torus, in the 



horizontal and vertical directions, respectively. The r.h.s. of Eq.(2.3) is 
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very similar to the multipolygon representation Eq.(2.1), modulo the sign 
(— ai) h<yT \— ct2) v ^ (— l) fi -( r ) t '( r ) ) which depends on whether the parity of the 
number of windings in the horizontal-vertical directions are even-even, or 
even-odd, or odd-even, or odd-odd. The value taken by this sign with differ- 
ent boundary conditions and different winding parities can be conveniently 
summarized in the following table. 





even-even 


even-odd 


odd-even 


odd-odd 


« = (+,+) 


+ 








« = (+,-) 


+ 


+ 




+ 


« = (-+) 


+ 




+ 


+ 


a =(-,-) 


+ 


+ 


+ 





More explicitly, we can write: 



Z° ++ ({J b }; A £ , L ) = Zl e ({J b }; A<, L ) - 

-Zl Q {{J b }- A lyL ) - ZleUM; All) - Z° _ ({J b }; A e>L ) 



(2.4) 



Z° + _({J b }; A i>L ) = Z° e _ e ({J b };A £)L ) + (2.5) 
+Z e _ o ({ J b }; A itL ) - Z o _ e ({J b }; Az, L ) + Z o _ o ({J b }; A LL ) , 

Z° + ({ J b }; A e>L ) = Z° e _ e ({J b }; All) - (2.6) 
-Zl_ Q {{ J b }; A ijL ) + Z o _ e ({J b }; V) + Z o _ o ({ J 6 }; A ltL ) , 

Z°_({J b }; A e>L ) = Z° e _ e ({J b };A e>L ) + (2.7) 
+Z e _ o ({ J b }; A e , L ) + Z o _ e ({J b }; A e , L ) - Z o _ o ({J b }; A £>L ) , 



where 



Zl e ({J b }; A LL ) = 2 £L [ n cosh(M)] Y^Zt L niI tanh (M) , (2-8) 

b£Bt, L ' 7er 6G7 

and the superscript (e — e) on the sum ^rcA* L indicates the constraint that 
T winds over A^l an even/even number of times in the horizontal/vertical 
direction, including the case that V does not wind over the torus; in other 
words, when we say "even number of windings" , we include the case of zero 
windings. Of course, the other partition functions, namely Z®_ ({ J b }; A^l), 
Z®_ e ({ J b }; A^l), Z®_ ({ J b }; A^l), are defined similarly, with the constraint 
that T winds up over the torus an even/even number of times in the hori- 
zontal/vertical direction replaced by the one that T winds up an even/odd, 



8 



odd/even, odd/odd number of times, respectively. By definition, the total 
partition function in Eq.(2.1) is 

Z°{{J b };A i>L ) = Z° e _ e {{J b };A tjL ) + 

+Zl ({J b };A^ L ) + Z o _ e ({J b };A e , L ) + Z° _ ({J b }; A,, L ) . (2.9) 



Alternatively, using Eqs.(2.4)-(2.7), we can also write 

Z°({J b }; A e>L ) = ^(Z°_({J b };A e>L ) + (2.10) 

+Z°_ + ({ J b }; A ijL ) + Z° + _({J b }; A t>L ) - Z° ++ ({ J b }; A E , L )) , 

which is the desired connection between the multipolygon and the Grass- 
mann representations. This relation is valid for all {J b } b& t3 lL and all inverse 
temperatures (3. On the other hand, if J b is independent of b, the identity 
Eq.(2.10) gives us a mean to compute the partition function in closed form, 
simply because the Grassmann integrals Z^(A tjL ) := Z^,({J b }; A^ L )| Jb=J are 
gaussian and translation invariant. In practice, the computation proceeds as 
follows: one first goes to Fourier space, thus block-diagonalizing the quadratic 
action appearing in the Grassmann integral; each block one is left with in- 
volves the degrees of freedom associated with the Fourier modes k and — k, 
with k e V a , and 

V a = \ k= ^-j(r+^^), T (n+^^)j : r = 0, . . . ,£-1; n = 0, . . . , L-l 

(2.11) 

The computation of the Grassmann integral of the variables associated with 
each block is elementary and leads to a determinant or to a Pfaffian, de- 
pending on whether k differs from — k or not. The result is the following. If 
a ^ (+, +), then, defining t = tanh(/3J) and S t ($) ■= S({J b }; $)!,,, 



Z° a (A. 



I V $ e s t m = "Q [(l+t 2 ) 2 -2t(l-t 2 )(cosA; 1 +cosA; 2 )] 1 / 2 . 

(2.12) 



(2cosh 2 (/3J)) tt 
If, on the contrary, a = (+,+), 

„%L = / V •«*<•) = (2 - (1 + tf)(2 - (1 - If) 
(2 cosh 2 ((3 J)) J 



Yl [(l+t 2 ) 2 -2t(l-t 2 )(cosk 1 + cosk 2 )} 1/2 , (2.13) 



kex>++: 

k^O,(7r,7r) 



which is positive for < f3 c , negative for (3 > (3 C and vanishes at the critical 
point (the critical point /3 C is defined by the condition that t = y2 — 1). 



9 



The difference in the results obtained for a ^ (+, +) or a = (+, +) is due 
to the fact that in the first case all the modes k £ D a can be grouped into 
pairs (k, — k) and, correspondingly, the evaluation of the gaussian Grassmann 
integral reduces purely to a product over determinants (each determinant 
being the integral over the variables of the modes k and — k). In the second 
case, all modes but two can be grouped into pairs, the two special momenta 
being k = and k = (vr,7r); therefore the result is equal to the product of 
the determinants associated with the paired momenta times the two Pfaffians 
coming from the modes and (tt, 7r), which give the factor (2 — (1 + 1) 2 )(2 — 



[1 - tf) in the r.h.s. of Eq.(2.13 



In conclusion, evaluating Eq.(2.10) at (3 C in the translation invariant case 



and using the fact that the (+, +) Grassmann partition function vanishes at 
criticality, we find: 



Z°(A e 



0=Pc 



1 



{Z°__(A e , L ) + Z°_ + (A e , L ) + Z° + _(A^ L )) 



Z° + _(A t 



zU^L)\l(i + z ° +iAe ' L 



+ 



Z°__(A £<L ) Z°_(A 



(2.14) 



with 



Z°-(^,l) 

z° + (^ L : 

ZUHl] 



0=Pc 



0=Pc 



2) iL [ ) (4 - 2 cos h - 2 cos k 2 ) 



1/2 



keT>- 



(V2) iL \ I (4-2COSA;! -2cosA; 2 ) 1 / 2 (2.15) 



keV- 



{V2) eL \ l (4-2COSA;! -2cosfc 2 ) 1 



/2 



Taking the logarithm at both sides of Eq.(2.14), dividing by the volume, and 
taking the infinite volume limit, we get the bulk term f^, which is given by 
Onsager's formula: 



lim logZ°(A Ai ; 

,L->oo IL 



(2.16) 



log 2 



dk 

(2^ 



log(4 — 2 cos k\ — 2 cos k 2 ) 



Using the notation of Eq.(1.4), we write the first finite volume correction to 
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the critical pressure in the form with 



— = lim lim 

6 i— >oo L— >oo 



V2L 



log(4 - 2 cos A* - 2 cos k 2 ) - t (/«, - - log 2) 



ke£>- 



+ 



limlim « log ri + + |f#4) l. (2.i7) 



provided this limit exists and is finite. In Appendix A.l we show that the 
limit in the second line is equal to zero. On the contrary, the one in the first 
line is non trivial and can be explicitly computed as follows. Taking the limit 
L — > oo first, we can rewrite the first line as 



lim f y 



e/ ±? y-fn+f dkl r dk 



J Sn~j 



2vr 



2tt 



where £ n 



2tt , 



[log(4— 2 cos£ n — 2 cos A; 2 )— log(4— 2 cos k\— 2 cos fc 2 

(2.18) 



n 



|). The integral over fc 2 can be performed explicitly [18, 



Formula 4.224(9)], leading to 



lim f y 



*/2-l „t r E „ 



=0 ■'e. 



2tt 



£/2-l 

WW-7(*i)] = lim e 



n=0 



7 dtf 



[l{in)-l{in + k r )\ 

(2-19) 

with 7(fe) := cosh 1 {2 — cos fc). Expanding in Taylor series 7(£ n + k') around 
£ n up to second order included, and using the fact that j(k) is a C°° function 
on [0, 7r], we find that Eq.(2.19) can be rewritten as 



lim 

l— >oo 



P 2 V 2 ' 1 



n=0 



(2.20) 



7T 

12 



7T 



^/(Jfe) = — ( 7 '(0 + ) - 7'(7r)) 



12 



7T 

12 



C7T 



which corresponds to c = 1/2, as desired. 



3 The interacting case 

We now attack the problem of computing the first non trivial finite volume 
correction to the critical pressure in the interacting, A ^ 0, case. We make 
use of the results and methods of [H], which we refer to for the proof of 
numerous relations used in the following. As in the previous section, we drop 
the dependence on in the symbol used for the partition function, and we 
assume that £ and L are even. 
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3.1 Grassmann representation of the interacting par- 
tition function 



The interacting partition function Eq. (1.2) can be written in a form analogous 



to Eq.(2.10), for all temperatures (3: 



Z(A e>L ) = i(2L_(A />2 + Z.+(A t , L ) + Z + -(A e ,L) - Z ++ (A e , L )) , (3.1) 
with Z a {h-t,L) given by (see P~H Proposition 1]): 

Z a {A t , L ) = C e>L j P$ e s ^ +v ^ , (3.2) 

where: 

• Cii, is a normalization constant, defined as 

C iiL = (2cosh 2 (pj)) iL e v ^ H cosh 2 (^(x - y)) (3.3) 

{x,y} 

with V^l{X) an analytic function of A, defined as (using the notation 
of [2], see the proof of [HJ Proposition 1] and, in particular, [HI 
Eq.(2.29)]) 

V^{\) = 21L £ ||^nC(M;7)- (3.4) 
supprgfco 

Here b is an arbitrary n.n. bond of Ail, supp T = LL 6 r7 and | supp T| 



is the number of bonds in supp T. The sum in Eq.(3.4) is independent 
of bo, by translation invariance; the activity £(0,0; 7) (defined in p~U 
Eq.(2.18)]) is a translation invariant exponentially decaying function, 
satisfying the bound (see [TU Eq.(2.28)]) 

|C(0,0;7)I<^, v = 4e^^(^) 1/Mo . (3.5) 

If we define E xA = H x H x+a&1 and £ Xj2 = V x V x+0 & 3 , then V($) is a 
polynomial in {E x j} x ~^ L , which can be expressed as 

n 

V($) = EE E W i {7L 1 ,...,x n )l[E Xitji (3.6) 

n>l ji,...j„xi,....x„ i=l 

12 



where j = (ji, ■ ■ ■ ,j n ) and the kernel Wjfai, . . . , x n ) is translation in- 
variant and satisfies the following decay Tjound: 



|H^(x 1? . . . ,x„)| < C7 n (/3|A|) max{1 ' cn} e- K5(xi ' 



•,X n ) 



(3.7) 



for suitable constants C, c, k > depending only on M (the range of 
the interaction); here 5(x 1; . . . , x n ) is the tree distance of the set X = 
{xi, . . . , x n }, that is the length of the shortest tree graph composed of 
bonds in Bi l which connects all the elements of X. 



The Grassmann integral Eq.(3.2) and the correlation functions of the $ field 
induced by the "measure" X>$ e St($)+v($) h ave b een studied in great detail 
in [Hj. In particular, part of the main result of [H] can be reinterpreted by 
saying that there is a critical temperature /3 C (A) such that, if we fix t = t c = 
tanh(/3 c (A) J), then the two-point function of the $ field decays polynomially 
(like distance -1 ) at large distances: if a. 7^ (+, +) and we perform the unitary 
change of variables from $ to the critical modes ip, \ defined as 



(3. 





A 




/#x\ 




/ e<i 


e l 4 1 —A 






= [/ 


#x 






e 4 1 % 


Xx, 


+• 


^x 




-e 4 


— e 4 1 — j 


\Xx, 


-/ 








■ 7T 


— e 4 ^ 1 z / 



then, asymptotically for large distances, 



Z(A) 



5 tc (<K)+V($) 



7Tt r 



(yi - zi) + - £2) ' 

(3.9) 

Here the symbol "~" means "up to faster decaying terms as |x — y| — > 00" 
and Z(X) is the analytic function appearing in [T^ Theorem 1.1]. An infinite 
volume limit £, L — > 00, performed while keeping the sites x, y fixed, is 



implicit in Eq.(3.9). Eq.(3.9) can be read by saying that the asymptotic 
behavior of (V'x.wVv^t is the same as that of (V'x.wV'yjW')^) z> where t° = 
\/2 — 1, Z = w^yfs and (-)^ t o is the average with respect to a properly 
renormalized massless gaussian integration: 



(3.10) 



A way of computing Z, similar to but slightly different from the one pro- 
posed in [2], will be described below. The integration T)Qe ZSt ° is the 
right reference measure, around which to perform the perturbation analysis 
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of Eg. (3.2) at t — t c , along the lines of [14] . 



However, before we start describing the renormalization group computa- 
tion of Z a (Ai^) at t = t c , let us discuss how to deal with the sign problem 
in Eq.(3.1): it is in fact apparent that the expression in the r.h.s. involves 
a difference between Grassmann partition functions, which may in principle 
produce dangerous cancellations between the different terms. This has to be 
contrasted with the computation in the A = where, as we saw above, at 
criticality Z+ + (Ai t L) = and the other three terms are all positive. This 
may even be true at A ^ and t = t c , but we do not know how to prove it. 
Nevertheless, we can prove an a priori partition function inequality, which 
can be thought of as a weak version of this claim, and is actually enough 
to the purpose of computing the pressure up to the first non trivial finite 
volume correction. This is discussed in the next subsection. 



3.2 A partition function inequality 



As we said above, our goal is to compute Eq.(3.1) at criticality, by us- 



ing the representation Eq.(3.2) and a renormalization group analysis for 
f T?$ e 5 *(*) +v (*), along the lines of [H]. However, as we will see below, 
our renormalization group computation of J P$ e St($)+v(#) on \j wor k s f or 
a. 7^ (+, +), due to the possible vanishing of Z ++ (A^l) at t — t c ; the problem 
is that at the unperturbed level Z+ + (Ae t L) = at criticality and, therefore, 
in order to perturbatively compute Z ++ (Ai tL ), we do not know where to per- 
turb around. Luckily enough, in order to compute the pressure up to the 
first non trivial finite volume correction, we do not really need to prove that 
Z ++ (Ai t i) is zero or much smaller than the other three partition functions: 
a weaker statement, summarized in the following Lemma, is actually enough 
for our purposes. 



Lemma 1. Under the stated assumptions on the potential v(x) in Eq.(l.l), 



for A > and all inverse temperatures (3 > 0, the following inequalities hold: 

i< Z ( A ^ <! f 3 in 

3 " 2_(A^) + Z^{A ljL ) + Z + _(A itL ) " ' 1 • ) 

Z_+(A £iL ) + Z + _(A, iL ) >0. (3.12) 

The proof of this lemma is postponed to Section |4| In order to use this 
result, we have to combine it with the information that 

2L_(JW) > , (3.13) 
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which is proved in the next subsection for A small enough. Putting Eqs. (3.11)- 



(3.13) together we get: 



\Z„ (Af L ) < 2(A«) < (A tl ) (l + fe^i + 
3 V Z — (A^lJ Z — (A£,lJ 

(3.14) 

Let us now put ourselves at criticality, t = t c . In Appendix A. 2 we prove 
that 



£ / |Z_ + (A,, L )| |Z + _(A <|L ) 
hm hm - log 1 + ^ — { + — — t-t^ 
^ool^coL V Z (A^i) Z (A^ iL ) 







(3.15) 



and, therefore, using the notation of Eq.(1.4) 



/«,= lim - 1 -logZ(A £iL )= hm -^-logZ__(A, iL ) , (3.16) 

£,Z,-»oo IcL £,L-*oc tL 



— = lim lim 

6 t— >oo L— S>00 



logZ(A £iL )-£ 2 /o 



lim lim 



logZ__(A,, L )-£ 2 /o 



(3.17) 

which reduces the computation of the central charge to the evaluation of the 

finite volume corrections to Z (A^l). These are computed in the next two 

subsections. 



3.3 Renormalization group analysis of Z (A^) 

We start from Eq.(3.2) with a. = (— , — ) and t = t c , and we rewrite it as: 



Z—(A e j,) = Ce iL / P$e 



ZS t o(*)+V($) 



(3.18) 



where 



V($) = v($) + S tc {$) - ZSiP c {$) . (3.19) 

In the following, t c and Z will be constructed in such a way that, asymptot- 
ically for |x| — > oo, 



(V'O,a)'0x, £ j')t - ('0O,wV'x,w')t° 



1 1 



5, 



Z ltt° c Xi + ilx!X2 



(3.20) 



that is the theory is critical (i.e. massless) and Znt c is the dressed wave 



function renormalization, see the discussion after Eq.(3.9). Equivalently, the 



t c and Z will be chosen in such a way that the flow of the running coupling 

constants defined in the iterative construction of Z (A^ L ) described below 

remain bounded at all scales; we refer to the following for the definition of 
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running coupling constants and for the study of their flow; see in particular 
the neighborood of Eqs.( 3.74~| )) and (3.76 for a discussion about how to fix 



t c , Z in terms of an implicit function theorem. 

It is convenient to multiply and divide Eq.(3.18) by the proper normal- 
ization, Z°__(A e>L ) = J V$e ZS ^\ thus finding 

Z— (V) = C ltL Z°__(A e , L ) J F(d$) e y W , (3.21) 

where P(d$) = [Z°_(A e , L )] e ZS ^ m is the normalized gaussian refer- 
ence measure. The normalization constants C^l and Z^_(A^l) are explicit 
and, therefore, the only non trivial part to deal with is 

S__(A^):= J P(d$)e V ^ . (3.22) 

In order to compute this integral we proceed by following essentially the same 
strategy of [2], modulo a few small modifications described below. As a first 



step, we pass to the critical modes, already introduced in Eq.(3.8); in order 
to fix the normalizations as in [TJ] we also rescale the variables as (see [T4"l 
Eq.(2.38)]) ip u — iu^irt® ip u , Xu> ~~ —ioJ^/irt® Xu» and next we perform the 
following linear change of variables (see [T4"l Eq.(2.52)])) 

Xk->Xk + ^ 1 (k)g(k)4, (3.23) 
where, if cr x (k) = cos k x + cos k 2 + 2 v ff~ 1 , 

C (k) = (~ l Sin kl + Sin k2 i(Tx ^ ^ (3 24) 

x I —ia x (k) — « sin /ci — sin A^y ' 

„ ,, v _ / — z sin ki — sin fc 2 « cos fci — z cos A; 2 \ 
L ' z cos fci + z cos &2 — z sin fci + sin k 2 J 

After these changes of variables we get the analogue of [2J Eq.(2.53)], namely 
H__(A,, L ) = / P{d^)P{d X )e V ^ , (3.26) 



where: 

• If, for k e P , we define 
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then the Grassmann gaussian integrations P(dip), P{dx) can be written 

as 

p(#) : =j-[ n n ex ^{~i^j e ^ T k^(k)^ k } 

:= [ It II ^ k '-] ex P { - ^Zl ? *-k^x( k )Xk} . 



kez> oj=± ke£> 

where, letting cr^(k) := cos/ci + cos/^ — 2, 

C,(k) = ^-^in^ + Bin^ ^(k) \ _ g(k)c - 1(k)g(k) 

^ y — «<j^(k) — i sin A/l — sm k^) x 

(3.27) 

and Af^,Af x two normalizations, such that / P(dip) = J P(dx) = 1. 

V(/0) x) is the rewriting of V($) in terms of the new variables. It is 
easy to check th at it s kernels satisfy the same decay estimates as those 
of V($), see Eq.Q- 



The propagators of the if) and x fields are, respectively, 

kex> 

^ 1 , / - Pi 



kex>_ 



A simple explicit computation shows that asymptotically as |x — y| — > oo the 

propagator of the ib field behaves like 4 7 v> as it should, while 

the one of the x field is exponentially decaying. Therefore, we can integrate 
out the x field, by proceeding as described in detail in [HI Section 3.1], and 
we get the analogue of [131 Eq.(3.1)]: 

H__(A^) = e L£E ° [ P(#)e v(0) W (3.28) 



where Eq is an analytic function of A, vanishing at A = and weakly de- 
pending on the volume; namely, if we denote by Eq 3 its infinite volume limit, 
then \E — E^°\ < C|A|e _rf for two suitable constants C, c > 0; see below for 
a more detailed discussion of the finite volume corrections to the thermody- 
namic quantities of the theory. 

Multiscale analysis. We now need to integrate out the if) field. However, 
since the if) field is massless, we cannot do it trivially in one step. A convenient 
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way to proceed is by an iterative procedure, described in detail in [TU Section 
3.2]. We define a sequence of geometrically decreasing momentum scales 2 h , 
with h = 1,0, .. . Correspondingly we define a sequence of cutoff functions 
//j,(k) in the following way. Let x{t) be a smooth compact support function 
that is 1 for t < 1 and for t > 2. We choose /o(k) = 1 — x(|k|) and 
/^(k') = x(2" /l |k|) - x(2- ft+1 |k|) V7i < 0, so that f h for h < is non zero 
only if 2 h ~ l < |k| < 2 h+1 , and 

l = ^A(k). (3.29) 

h<0 



The resolution of the identity Eq.(3.29) induces a rewriting of the propa- 



gator of ip =: ^>(-°) as a sum of propagators concentrated on smaller and 
smaller momentum scales and an iterative procedure to compute Z. At each 
step we decompose the propagator into a sum of two propagators, the first 
approximately supported on momenta ~ 2 h (i.e. with a Fourier transform 
proportional to /^(k)), h < 0, the second approximately supported on mo- 
menta smaller than 2 h . Correspondingly we rewrite the Grassmann field as 
a sum of two independent fields: ip(- h ^ = ipW + and we integrate out 

the field ip( h > in the same way as we did for \- The result is that, for any 
h < 0, we can rewrite 

H__(A, L ) = e UE * / P(#(^))e vW «' ( -' i) ) (3.30) 



where E^, are defined recursively, and P(dip ( -- h ^ > ) is the gaussian integra- 
tion with propagator given in momentum space by Ylh'<h 9 CO) with 

yW(k) = yA(k)[C^(k)]- 1 . (3.31) 

Note that the direct space counterpart of g^ h > decays to zero faster than any 
power: 

\g (h) (x)\ < 7 j? c/ \ i \ , Vp > . (3.32) 

|y v ;| ~ (1 + 2 h \S(x)\)p v ; 

where 

5(x)=(-sin(— ),-sm( — )). 

The outcome of the iterative construction is that, in particular, the effective 
potential V^ h ' has the following structure: 

2n 

v (h) w = ee ^ ( ;L(x)[n^^] > (3.33) 

n>l iii,x i=l 
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where x = (xi, . . . ,x 2n ), u = (ui, . . . ,u 2n ) and 

\- 9 E l^ 2 ( n;l(xi,...,x 2n )|<C'' 1 2 h ( 2 -")|A| max ^ c "> (3.34) 



Li 

Xl,...,X 2 „ 



for suitable C, c > 0. For future reference, let us also rewrite Eq.(3.33) in 
momentum space: 



2/i 



vW W = E?zirT E ^(ki,...,W 1 )[n^ lWl ]«5(k 1 + ---+k aB ) 

n>l ^ > ki,...,k 2 „ i=l 

(3.35) 

where <5(k) is a Kronecker's delta, periodic over the torus IR 2 /27rZ 2 . 

The iteration continues until the scale h* := [log 2 (tt/^)] is reached. At 
that point, the left-over propagator, g(- h *"> is massive on the "right scale" 
(i.e. on the very same scale 2 h *), so that the associated degrees of freedom 
can be integrated in one step. The result is the desired partition function. 



Localization and renormalization. In order to inductively prove Eq.(3.30) 
we write 

y W ty) = CV {h) (VO + nV {h) , (3.36) 

where £V {h) (ip) and T^V^(^) are the so-called local and irrelevant part of 
the effective potential, defined in the next few formulas. We use a definition 
of localization operator on the lattice at finite volume, analogous to the one 
used in [3] where the finite volume effects of the renormalization procedure 
are discussed in great detail. To be fair, the definitions below are a bit more 



complicated than those in j3] (see the slightly cumbersome definitions (3.38) 



( 3.39[ )), the reason being that we want to make sure that the relative finite size 



errors induced by the localization procedure are of the order 0[l 4 ) rather 



than Q(£ 2 ), wh ich would not be enough to our purposes; see Eqs.(3.42) 



(3.44) and (3.48) below. Let us now come back to the definition of £V^(ip) 

■(h) 



and TZV^(ip). If we think of the kernel W^ o.fw w \ as a 2 x 2 matrix with 



matrix indices oji,U2, we let 

CVWy,) = -^E^[^2 (/l) (k)]^-k (3.37) 

k 

1 ' ki,k 2 ,k 3 
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where, setting k w / = (7y|,7/^,, 



and 



77,77' =± 



9, 



Wf } (k w )(l + ^ 



sin fci , sin &2 



+ 7?" 



sin(7r/£) sin(7r/L) 



-^ 2 W (3k w )(l + ^ 



sin fci , sin hi 



+ rf — 



sin(37r/£) sin(37r/L) 



£W^( kl ,k 2 ,k 3 ) = £ 



-W^(V / k / k / 



(3.38) 



(3.39) 



Note that in the limit L,£ — > oo, the action of the localization operator 
reduces to: 



£Wf } (k) = W 2 W (0) + (sinfei^ + sinfc 2 a fc2 )W 2 w (0 



r{h) 



and, similarly, 

£^(^,^^3) = # 4 ( fJ(0,0,0) . 

In other words, at finite L and £, CW 2 (k.) and £1^4.^; (kj, k 2 , k 3 ) have to be 
understood as finite volume, lattice versions of the Taylor expansion up to 
order or 1, respectively. The coefficients 9/8 and —1/8 have been fixed in 
such a way that the difference between the finite and infinite volume local- 
ization operators goes to zero as [max{L, £}]~ 4 ; if needed, we could further 
modify the definition of localization in such a way that the operator remains 
the same in the infinite volume limit, but the difference with its finite volume 
counterpart goes to zero faster than any power of max{L,£}. Note also that 
C has the nice feature of being a projection operator: C 2 = C 



The direct-space counterparts of Eqs.( 3.38 )-( 3.39) read as follows: 



(<*) 
y 



(3.40) 



£ ^- h W\^ y) \Gy,^ h) + E *(y» x )<^ 



x,y 
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and 



£ e ^( x i> x 2,x 3 ,x 4 )n^i = 

Xl,X2,X 3 ,X 4 j = i 
UJl ,UJ2 ,^3 ,^4 

4 

E <2(*i, * 2 , x 3> x 4) II [gwSS 



(3.41) 



Xi,X2,X3,X4 
,tx>2 ,^3 ,W4 



i=l 



where G y , x and dj(y,x) are translation invariant and 

„ 9 /7TXix /7rx 2 x 1 /37rxis /37ra;2N 
G x , = - cos (— ) cos ( — ) - - cos ( — ) cos ( — ) 



(3.42) 



9sin(gf) vrx 2 , 
di(*, 0) = - , \;/ cos (— ) 



sin 



(!) 



lsin(^i) 
(¥T 



sin 



cos (^) , (3.43) 



d 2 (x,0) 



9 /7 rxi N sin(^) 1 ^ gin (^) 
-cos( — -) cos(— — )■ 



sin 



(!) 



sm 



(?) 



(3.44) 



and di is the symmetric discrete derivative w.r.t. Xi, i.e., d\f(x\,X2) = 
2if{ Xl + l' x 2) — f(xi — 1,^2)], and similarly for <9 2 . A few comments are 
now in order. 

1) The action of C in direct space can be seen as an action on the fields, as 
indicated by Eqs. ( 3.40 )-( 3.41). The action of 1Z = 1 — C can be interpreted 



in the same way, too: it takes the form 



(3.45) 



* E & n v>g2 = E < & ^giU^g iict^gs + 



i=l 



^ ^ xi ,x 4 S^x 4 ,wi ^x 2 ,x 4 ;oj 2 ^x 3 rx 4 ,oj 4 ^ ^ xi ,x 4 V y X4 ,wi ^ x 2 ,x 4 Y x 4 ,oj 2 x 3 ,x 4 ;o; 3 r : 



X4,0J4 



where 



r£g = 4% h) - G y , x ^ fe) - E *(y. x W 



£>(</>) = -/,(</») _ £ ./,(</») 



(3.46) 
(3.47) 



2) The functions G yjX and di(y, x) are antiperiodic over A^x in both their 
arguments. This implies that the fields T y ^2 and -D y ~x^L are antiperiodic 



in y and periodic in x. Therefore, the summands in Eqs.(3.40)-(3.41 ) are 
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periodic both in y and x. 

3) For fixed x, asymptotically as L ^> £ ^> |x|, we can rewrite 

G x ,o = l + 0((|x|/£) 4 ) , ^(x,0)=x,[i + O((|x|/£) 4 )] . (3.48) 

which is useful for bounding dimensionally the effect of the finite volume on 
the localization procedure. 

4) The action of 7Z produces a dimensional gain on the Grassmann monomi- 
als which it acts on non-trivially. This can be seen as follows. As discussed in 
item (1), the action of 1Z can be thought of as a replacement of a field ijj y ~^ ei- 
ther by or by D y -£\ with x, y G A iiL (recall that the box is topologically 
a torus); since both x and y are represented by infinitely many equivalent 
images on Z 2 , it is always possible to pick two such images, to be called x, y, 
so that their euclidean distance on Z 2 is the same as the distance between x 
and y on the torus A^ L . We have D y jH ] = where 



D 



(<h) 



Gy^ipsf^ can be conveniently written as: 



D 



(<h) 
y,x;o; 



(y 



(<h) 
x' 



x'=x+s(y— x) 



+ (1 - G y ^ 



(<h) 



(3.49) 



where (1 — G y , x ) is of the order 



|y - x| 4 ^ 



Note that when we plug this 
formula in the last two lines of Eq.(3.45), the factors (y — x) or (1 — G y ^) 
multiply the kernel W^ h \ whose value is obtained by integrating the degrees 
of freedom on scales strictly larger than h, i.e., it is a sum over Feynman 
diagrams with propagators g^'fa — y«) all of scale hi > h. Consider the 
contribution from the interpolated term first. When we decompose the factor 
(y — x) along the diagrams contributing to W^ h \ some of the propagators 
g( hi \~x.i — y~i) are multiplied by (x» — y,), which is dimensionally equivalent 
to a factor 2~ hi (see Eq.(3.32)). On the other hand, the derivative d x > inside 
the interpolation integral, when acting on ip^ h \ is dimensionally equivalent 
to a factor 2 h , simply because is a field that is smooth on scale 2~ h . 



Similarly, the contribution from the second term in Eq.(3.49) can be bounded 

4 and then by decomposing 



xi 4 r 



by first replacing (1 — G y ,x) by (const.) |y 

the factor |y — x| 4 along the diagrams contributing to HQ , so obtaining 
2-4h,£-4 _ 2&{h*-hi) f or some > ^ j n conclusion, the action of 1Z on the 

four-legged kernels is dimensionally equiv alent to a dimensional gain 2 h ~ hi , 
with hi > h, if we pick the first term in Eq.(3.49), or to 2 4( ^ h *~ h ^ if we pick the 



second term, which is due to the finite volume corrections to the localization 
procedure. A similar discussion applies to the action of 71 on the two-legged 
kernels, in which case 1Z is equivalent to eith er 2 2 ( h ~ hi \ with hi > h, if we 
pick the analogue of the first term in Eq.(3.49), or to 2^ h *~ hz \ if we pick the 
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analogue of the second term. 

5) A key fact which makes the theory at hand treatable (and asymptotically 
free) is that the action of C on the quartic kernels is zero "by the Pauli 
principle", i.e., simply by the Grassmann rule ifj^w = ®- m f ac t ; note that in 



the second line of Eq.(3.41) at least two of the four u> indices must be equal 



among each other. Therefore the expression in the second line of Eq.(3.37) is 



identically zero. This property can be diagramatically interpreted by saying 
that the fermionic nature of the theory automatically renormalizes the four- 
field interaction, which is dimensionally marginal (see below) but effectively 
irrelevant thanks to the cancellation that we just mentioned. 
6) By the previous comment, the only non vanishing contribution to the 
local part of the effective action is the one in the first line of Eq.(3.37). By 



its very definition, it is apparent that £W" 2 is invariant under the discrete 
symmetries of the theory (reflections, discrete rotations, etc, see [TH Section 
II. D] for a thorough discussion of this point). Therefore, its most general 
form is (see [H (2.68)]) 



£Wf } (k) 



Zh(i sin k\ — sin k 2 ) 



i2 h v h 
Zh(i sin ki + sin k 2 



(3.50) 



for two real constants Zh^h, which will be called the running coupling con- 
stants of our theory. The initial values of these constants Zq,i> q , are induced 
by the choices of t c , Z in Eq.(3.19), as well as by the effects of the change of 



variables described after Eq.(3.22) and of the integration of the massive fields 



in Eq.(3.28). It is straightforward (if lengthy) to keep track of this series of 



transformations and to check that zp, vp are analytically invertible functions 
of t c , Z in a neighborood of t c — t°, Z — 1, A = 0: 



>o, fo) = {F (X, t c , Z), iV (A, t c , Z)) & 



Z = l + ((\,z ,vp) 
-. = t° c + r(A, z , u ) 



(3.51) 



where F , N , (, t are analytic functions of their arguments. It is straightfor- 
ward (if lengthy) to check that, if |A| + \uq\ + \zq\ < £o, 



dz 
drjX, go, up) 
dz 



a + O(e ) 
O(e ) , 



d((\, z , up) 



dvp 

dr(X, z , up) 
dvp 



O{e ) , (3.52) 



b + O(e ), (3.53) 



for two non vanishing constants a and b. This property will play an important 
role in the following, in the choice of the parameter Z. 
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Tree expansion. Going back to the inductive proof of Eq.(3.30), we note 



that representation is valid at the first step, see Eq.(3.28). Assuming the 



representation to be valid at scale h, let us show that the same structure is 



preserved at the following step. By using the rewriting Eq.(3.36) and the 



addition principle (see e.g. [13l Section 4]), we rewrite Eq.(3.30) as 



-UE h 



P(dip 



(<h-l)\ 



P(dip 



(3.54) 



where P{dip is the gaussian integration with propag ator^( h )(k),see Eq.(|3.3l[) 



By integrating out the degrees of freedom on scale h, which are massive with 
mass of the order of 2 h , we get: 



,Ll{E h +e h ) 



(<h-l) n V(' 1 - 1 )(V ( -' l - 1) ) 



P(d^ Sh - l) )e 



(3.55) 



Lleh + V^W) = log J P(# 



{h) \ £VW (ip+ipW )+KV (h > 0/M-V> (h) ) 



which proves the inductive hypothesis Eq.( 3.30[ ), provided that E h _i is fixed 
as 

£? h _i = E h + e h = E + e i • ( 3 - 56 ) 

h<j<0 

Note that the above procedure allows us to write the running coupling con- 
stants Zh, Vh with h < 0, in terms of Zk, Vk with h < k < 0: 



z h-i — z h + Ph 



v h _ x = 2v h + B 



h ■ 



(3.57) 



where df = 8f((zh, fh), ■ ■ ■ , {zo, ^o)) is the so-called Beta function. 

By the very definition of truncated expectation (see e.g. [TU (2.57)]), the 



second line of Eq.(3.55) can be rewritten as 



Lle h + V {h - l \*P) = 



(3.58) 

■CV {h \ip + -)+KV {h \i) + -)) 



n>l 



n times 



where is the truncated expectation associated with the gaussian integra- 
tion P(dip^). Iterating this relation, we are led to a tree expansion for the 
effective potential, as described in detail in e.g. [TU Sections III. A- III.D]. 
The resulting trees are defined in a way that is very similar to those in [T4"l 
Sections III.D], with the following minor differences (we refer for comparison 
to the description of trees in items 1 to 5 in [HJ Section III.D], to be called 
GGM.1-GGM.5, as well as to the discussion following that item list): 
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1) The trees we consider have only normal endpoints: m = in the notation 
of item GGM.l. 

2) The ultraviolet scale N (in the notation of GGM.2) is replaced by and we 
denote by 7n the set of labeled trees with root on scale h and n endpoints. 

3) The first four lines of GGM.5 are replaced by "With each normal end- 
point v on scale 2 we associate a factor V(tp^-°\ x); here x '■= i'^ should be 
thought of as the field on scale 1. With the endpoints on scale h v < 1 we as- 
sociate either a factor 2 h u h ^ £) k)W ^V'k.wV'-k.-w or z h jj £] ka) ^ w {i sin k x - 
ui sin fc 2 )V'-k,a;, in which cases we shall refer to the endpoint as being of type 
v or z, respectively." 

In terms of the definitions of trees, the effective potential can be 
written, in analogy with [HI (3.44)] as 

L£e h + V {h) (ip { - h) ) = J2J2 vW ( r ' ^) ' ( 3 ' 59 ) 

where, if v is the first vertex of r, if Ti, . . . , r s (s = s„ ) are the subtrees 
of r with root t>o, and if S^ +1 is the truncated expectation associated to the 
propagator g( h \k), 

(3.60) 

and F (?t+1) (r 4 ,^(^ +1 )): 

• is equal to TZV^ h+1 \ri, ip(- h+1 ^)) if n is non trivial; 

• is equal to 2 h+l Uh+ijj Xlk w ^V'k.wV'-k.-u; ^ r * * s trivial, h < N and the 
endpoint of n is of type v\ 

• is equal to Zh+ijj, ^^4^^ sin fci — u> sin k 2 )-ip-^ UJ if t* is trivial, h < N 
and the endpoint of r» is of type z\ 

• is equal to V(^- N+1 ') if Tj is trivial and h = 0. 

The values of the trees can be estimated dimensionally as described in de- 
tail in [HI Sections III.D.1-III.D.2]. In a notation analogous to the one 
introduced in those sections, we get the analogue of (T4, (3.35)], namely, if 
Wt,p,t,/8 is the kernel of the renormalized effective potential labelled by a tree 
t G Tn h \ a set of field labels P, a spanning tree T and a set of interpolation 
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parameters /3, 



re E E \w T prAO\ <c n [u ( C \x\r^ p ^ 

/3eB T x„ () 



v e.p. 



. 2 M2-|I^ I) J JJ 2^-V)(2-!l^l-.*0P« 



(3.61) 



d not e.p. 



with = 25| Pi; | i 2 + £|p„|,4- 



Remark. The bound Eq.(3.61) is essentially dimensional, i.e., it can be 



understood (modulo the combinatorics) by replacing all the propagators and 
the integrations involved in the definition of the kernel by their dimensional 
estimates, that is |p^(x)| < (const. )2 k and J dx\g^(x)\ < (const. )2~ k , see 



Eq.(3.32). To be fair, the dimensional estimate that we would obtain by 



these simple replacements would be similar to Eq.(3.61), but with z(P v ) 
replaced by zero. In order to justify the presence of the dimensional fac- 
tors n^riot e P 2~( hv ~ hv '' z ( Pv ' , one needs to take into account the action of 
TZ, see the discussion in item (4) above in the lines preceding and following 



Eq.(3.49). It is important to note that in order for the bound in Eq.(3.61) 
to be summable over the scale labels, it is not really necessary to have 
z(P v ) = 25\ Pv \2 + S\p v \,4'- a gaining factor z(P v ) = (1 + e)6\p v \ t2 + ^\p v \a 
with e G (0, 1) would make the job, too. For this reason, in order to regu- 
larize the kernels it is actually enough to use a portion of the gaining factors 
described in item (4) above; if desired, we can keep some of them on a side. 
E.g., concerning the factors 2^ h *~ hi ^ coming from the finite volume correc- 
tions to the definition of 71, we can use a fraction 2^ 1+£ ^ h *- h ^ with e G (0, 1) 
to regularize the kernels, and we can keep 2^- £ ^ h *- h ^ on a side. 



The bound Eq.(3.61 ) is valid provided that the running coupling constants 



rem ain o f order A, for all scales between h and 0. Under this assumption, 
Eq.(3.61 ) implies the analyticity of the kernels of and the decay bounds 



Eq.(3.34). 



An immediate corollary of the bound Eq.(3.61 ) is that contributions from 
trees r G 7n with a vertex v on scale h v = k > h admit an improved 



bound with respect to Eq.(3.34), with an extra dimensional factor 2 e ^ h k \ 
< 9 < 1, which can be thought of as a dimensional gain with respect to 



the "basic" dimensional bound in Eq.(3.34). This improved bound is usually 



referred to as the short memory property (i.e., long trees are exponentially 
suppressed); it is due to the f act th at the renormalized scaling dimensions 
d v = 2 — \\P V \ — z(P v ) in Eq.(3.61) are all < —1, and can be obtained by 
taking a fraction of the factors 2 (flv V)^ associated to the branches of the 
tree r on the path connecting the vertex on scale k to the one on scale h. 
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Under the same assumptions, the beta function itself, is analytic 
and dimensionally bounded by a constant independent of h. Moreover, the 
contributions to it from trees that have at least one node on scale k > h 
is dimensional bounded proportionally to 2 e ( h ~ k \ with < 9 < 1. It is 
remarkable that thanks to these bounds, the dynamical system induced by 
the beta function can be fully studied and shown to lead to a bounded flow 
of the running coupling constants. 

The flow of the running coupling constants. As announced above, the 
flow equations for the running coupling constants are 

= z h + f3 z h , = 2v h + dl , (3.62) 

where j3f = (3f((zh, u^), ■ ■ ■ , (z , z/ )) is an analytic function of its argument, 
with an analyticity domain bounded by: \zk\ + \v^\ < Eq, for all h < k < 0, 
where e is a suitable (small) positive constant. 

Note that both (3^ and f3% can be expressed as sums over trees with at least 
one endpoint on scale 2, the reason being that the local part of the trees with 
only endpoints of scale < 1 is zero by the support properties of the single-scale 
propagators that enter the definition of (3f. Therefore, by the short memory 
property, \/3 z h \, \0%\ < C e \\\2 eh , for 9 e (0, 1), uniformly in L,t The idea is to 
first solve the flow equations in the L, £ — > oo limit, by properly choosing the 
initial data in such a way that the sequence {(zh, Vh)}h<o remains bounded. 
Then we will use the same initial data as the L, £ = oo case in the finite 
volume equations and we will show that the resulting flow remains bounded 
and close to the infinite volume one, with explicit bounds on the error terms. 

Let us then consider the case L,£ = oo first. We denote by the 
corresponding beta function. We define 97tx,6> to be the space of sequences 
V. = (z,v) = {{z h ,Vh)}h<o such that \z h \ + \v h \ < K\X\2 9h , \/h < 0; we 
shall think of 9JIk,8 as a Banach space with norm || • \\g, where \\v\\e = 
su Pk<o(\ z h\ + Wh\)2~ 0k - Note that every exponentially decaying solution to 
the beta function equations (if any) can be looked for as a fixed point of the 
map T : WIr,6 —> ^-K,e defined by 

(Tz) h = -J2PP Z (v) , (TV). = -J2 2j ~ h ~^P l/ ^) ■ ( 3 - 63 ) 

3<h j<h 

The fact that, for K sufficiently large, T is a map from '!Mn;K,8 to itself is 
a simple consequence of the bound < Cg\\\2 eh . Moreover, if v,v' G 

yjlK,e, then using the short memory property: 

\PP*(v)-pp*(v!)\ < C' e \\\2^W-^\ < CM2^\\v-v[\\ e , (3.64) 

k>j 
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which implies that 



\\Tv-Tv\\ e < C's\X\\\v~v\\ e , (3.65) 

i.e. T is a contraction for |A| sufficiently small. Then the Banach fixed point 
theorem implies that T has a unique fixed point v* in WIk,b, which represents 
an exponentially decaying solution to the flow equations, with initial data of 
order A and given explicitly by the following expressions: 

z* = z* (\Z) := > "o = K(KZ) :=-Y,^^T U ^) ■ 

(3.66) 

By the previous construction, the functions Zq(X, an d u o(\ are analytic 
in their arguments in a neighborood of A = 0, Z = 1. 

We now consider the case of finite L, I and we pick the same initial datum 
(zq, Vq) for the flow equation as for the infinite volume case: 

z Q = z*(X, Z) , u = Z) . (3.67) 

Denoting by v — {(zh, Vh)}h*<h<o the sequence of running coupling constants 
generated at finite volume by this initial datum (using the notation above, the 
infinite volume counterpart of this sequence is denoted by v* = {(z^, f£)}/i<o), 
we get 

Based on this equation, and thinking of v as an infinite sequence (obtained 
e.g. by posing Vh = 0, Wh < h*), we can prove that, for any 9 < e < 1, 

\z h - z* h \ + \v h - u* h \ < C e , £ 2 eh ( 2 -p]' £ . (3.69) 



The proof is by induction in h. If h = the claim is obviously true, simply 
because the l.h.s. is zero. For h < 0, assuming the estimate to be valid for 



all the scales h < k < 0, we use Eq.(3.68), by rewriting the expressions in 
square brackets as 



W(v) - fi°'*m + w*m - /?rV)i . (3.70) 



Now, the second term can be bounded as in Eq.(3.64), |/3°°' # (w) — (3°°'*(v*)\ < 
C' e \\\2i e J2k>j \vk — v t\i so that by using the inductive assumption we get that 
the corresponding contribution, J2h<j<o l/^°' # (^) — 0j° (v*)\, is bounded 
from above by C'^\\\2 2h9 (2- h i~ x f'\ as desired. 
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The first term in Eq.(3.70) is due to the finite volume corrections. Re- 
member that, by using the tree construction explained above, Pf(v_) can be 



written as a sum of the form jj 



y y 



W, 



*< 



x nt J, where a is 



a suitable multi-index (collecting the indices P,T, f3 indicated in Eq. (|3.61 )) 
and W is periodic over A^l in all its coordinates. By construction, Wjjx) is 
a combination of propagators on the scales indexed by the tree labels, as well 



as of the functions G, di in Eqs.( 3.42 )-( 3.44) resulting from the action of 1Z 
on the nodes of the tree. Moreover W^{xi, . . . , x n J is translation invariant, 
so that we can fix one variable to 0, for instance Xi, and write: 



fSf (E) - /3~'*(E) 

= E[ E 

T,a x 2 ,...,x nQ eAf, i 



oo,#. 



(3.71) 



Wf (0,x 2 



E 



, x r 



X2,- 



where W"" '* is the infinite volume limit of Wj , which differs from the latter 
because of the replacement of the factors G and di by 1 and x iy respectively, 
and for the replacement of the propagators g^ k ' by their infinite volume limit. 
Now, the two sums YE 

' ^— 'X2,...,X rl 



X 2 , 



in Eq.(3.71) can be written as yi_ 2 
< £/A, Vi = 2,...,n' a 



, where V* v is over the set Ix, 
involves at least one coordinate outside the ball B, 



x + 

while 



£/4}. The easiest terms to bound are 



{x : |x| < 



E 



Wf (0,x 2 . 



+ |iy°°' # (o,x 2 



which are of the order |A|2 je (2 i /£) p for an arbitrary p > 0, simply because 
W* 



; contains a chain of propagators (each decaying faster than any power on 
a scale hi > j, see Eq.(3.32)) connecting with a coordinate x outside the 
ball Bi 4 . We are left with bounding 



E E (°> x 2 . . . , x n j - ^°°' # (o, 



(3.72) 



r,a x 2i . 



^/4 



where the differences in square brackets can be written as a sum of terms 
each of which involves either the difference between G and 1 (or, similarly, di 
and Xi), or the difference between a propagator and its infinite volume 
limit. Regarding the first class of terms, remember that the relative differ- 
ence between G, di and their infinite volume limit is bounded dimensi onall y 
by 2 4 ( h *~ hi \ see the discussion in item (4) above, in particular after Eq.(3.49). 
Part of this factor (at least a portion 2^ l+£ ^ h *^ h ^) is needed in order to renor- 



malize the kernels of the effective potential, see Remark following Eq.(3.61) 
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What we are left with is exactly a factor 2^ 3 s '( h * hi \ as commented at the 
end of that Remark. Summing these contributions over the scales j > h gives 



the desired bound Eq.(3.69). Finally, regarding the terms involving a differ- 



ence between the finite and infinite volume propagators, we observe that by 
the Poisson summation formula, the finite volume propagator (^(x), which 
is antiperiodic in its argument, can be written as a sum over images: 



nez 2 



1 



\ n ^g^(x 1 + n 1 £,x 2 + n 2 L) 



+ Sg 



CO i 



(3.73) 

where 5g^ h \x) is smaller than any power of (2~ h /£), namely, if |x| < £/2, 
\5g {h \x.)\ < C p 2 h (2- h /£y, \/p > 0. Therefore, the terms involving a differ- 
ence between a finite and an infinite volume propagator are smaller than any 
power in (2~ J /£), which is more than enough to the purpose of deriving the 



desired bound Eq.(3.69). This concludes the proof of that bound. 



In order to complete the discussion related to the choice of the initial 
data z , u , we are left with inverting the relation for Z, which is obtained 



by combining Eq.(3.67) with the infinite volume limit of Eq.(3.51), namely: 



l + ((\,z*(\,Z),v*(\,Z)) 



(3.74) 



The key ingredients to be used are the derivative estimates Eq.(3.52 ) together 



with the observation that the derivatives with respect to Z of the propagators 
obey to the same decay bound as the propagators themselves, so that 



84 
dZ 



-E 

i<o 



dZ 



0(A), 



(3.75) 



simply because dzl3°°' z (v*) = 0(X2^ e ') for some 9' < 9. Therefore we can 
apply the implicit function theorem to invert Eq.(3.74). This concludes the 



discussion about the choice of the initial data for the flow equation and, 



correspondingly, of the parameter Z in Eq.(3.18). Note that once Z is fixed 



the critical temperature t c is given by Eq.(3.51): 



t c = t° c +T(\,Z*(\,Z),V*(\,Z)) 



(3.76) 



3.4 Explicit computation of the bulk and finite volume 
corrections to the pressure 

In this section we compute the bulk and correction terms from the free en- 



ergy at the critical temperature, on the basis of Eqs.(3.16)-(3.17) and of the 
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construction of Z (A^) in the previous sections. Using Eq.(3.21 ), we write 

log (A ljL ) = log C e>L + log Z°_ _ (A tiL ) + log / e v ^ (3.77) 



and we compute separately the contributions from the three terms. 



The term logC^L- Using Eq.(3.3) and the properties of the potential v(x) 
spelled after Eq.(l.l), as well as the definition of V^l(A) in Eq.(3.4), we can 
write 



logC i>L = log(2cosh 2 (/3 c J)) + ^logcosh(^ A i;(x)) + 



x£2 



+ 2 £ 



suppTgbo 



IIc(0>0;t), 

7er 



(3.78) 



supp r| 

where /3 C is the interacting critical temperature (such that t c = tanh(/3 c J), 
with t c fixed as in Eq.(3.76)) and &o can be chosen aribitrarily, e.g., it can 
be fixed to be the bond connecting with ei. Note that the last term in 
the first line is independent of £, L, because u(x) has finite range, while the 
term in the second line differs from its infinite volume limit by exponentially 
small terms, which correspond to the contributions from multipolygons T 
that either wind up over the torus l or touch the complement of A^l on 
Z 2 (their exponential smallness follows from the decay bound Eq.(3.5)). In 
conclusion, the term log C^i contributes to the bulk term of the free energy 

(3.79) 



Eq.(p6b, 



i 



oo;l 



lim f-logC,, L = log(2cosh 2 (/3 c J)) + 



+ ^logcosh(^(x))+2 J2 



TCZ 2 : 
supp T3bo 



supp r 



IIC(0>0;7) 



7er 



but not to the finite volume correction Eq.(3.17). 

The term logZ°__(Ai l). By its very definition, see the lines preceding 
Eq.(3.21), we can write 



£L 



logZ°__(A e , L ) 



£L 



log / £>$e 



2\ogZ+ — log 



s t o(*) 



(3.80) 

where in the last identity we performed the Grassmann change of vari- 
ables $ -> VZ<$>. The last term, i log / £>$e 5t ? (<1>) , is the non interact- 
ing pressure evaluated at the critical point (up to an additive constant 
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jjr log Ce t L | A=0 p = p = \/2 + l), whose bulk and finite volume corrections have 
been evaluated in great detail in Section [2} Putting things together we get: 



/oo ;2 := lim — logZ°_(A, iL ) 



(3.81] 



, 1 

log 72TT + 2 



dk 



log(4 — 2 cos fci — 2 cos £2) 



and 



lim lim 



logZ°_(A £iL )-£ 2 /oo;2 



7T 

12 



(3.82) 



The term log J P(d$) e V( -*- ) . By the Renormalization Group analysis of 
J P(d&) e v ^ described above, we can write: 



tL 



log / P{0) e 



v(») 



(3.83) 



h*<h<0 



with £"0 and defined by Eqs.(3.28) and (3.55). Using the fact that + 
l^jfcl < C|A|2 and the tree expansion explained above, we can bound 



\e h \ < C\X\2 2h 22 h 



(3.84) 



uniformly in £,L, for all 9 G (0,1). Denoting by Eq 3 and e£° the infinite 
volume limits of Eo,eh, we can write the contribution to the bulk free energy 
under consideration as: 



/oo-s := lim log / P(dd>) e v ^ = JSg° + V < 

^ /i<0 



(3.85) 



which is an exponentially convergent series, whose sum is of order A. The 
finite volume correction of interest can then be written as 



lim lim 

l— >oo L—>oo 



log / P{d$) e v($) - £ 2 f ( 



oo;2 



(3.86) 



lim t 2 [(E - E™) + & - o - E 



h*<h<Q 



h<h* 



where E Q = lim^oo E and = linic-^oo e^. Now, (£?o — -Eo°) is expo- 
nentially small in £, as it follows from the fact that it can be written as 
a sum of terms that involve at least one difference between the finite and 
infinite volume propagator of the x field, which is exponentially small in l\ 
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therefore lim^ t>0 ^ 2 (£'o — E™) = 0. The last contribution to the r.h.s. of 



Eq.(|3.86|) is also easy to estimate; in fact, using Eq.(|3.84|) and the definition 

(3.87) 



of h* = Llog 2 (7r/£)J , we get: 

K\ < (const.)|A|2 2/l *2^* < (const.)|A|r 



2-g 



h<h* 



which implies that lim^oo i 2 



h<h 



0. We are left with 



lim I 2 



(3.88) 



h*<h<0 

Remember that (as well as eh and e^°) is a function of the whole sequence 



of coupling constants v (see the lines following Eq.(3.67) for a definition of 



v and of its infinite volume counterpart v*); we shall also indicate by v the 
L — > oo limit of v. We can then write Eq.(3.88) as 

}im£ 2 £ - + - e^*))] . (3.89) 



h*<h<0 



Thanks to the bound Eq.(3.69) on the difference between the finite and infi- 



nite volume running coupling constants, we immediately see that the second 
term in square brackets can be bounded as: 



\e?(v) 



< (const.)|A|2 (2+ 



i)h 



2- h \3 



(3.90) 



Picking 9, e > small enough, we get 



h*<h<0 



e?(v)-eW 



< 



(const.) |A|2^ 



3-e 



0. 



< (const.) | A | 



(3-91) 

We are left with the first 



which implies lim^oo £ 2 (e£° (v 
term in square brackets in Eq.(|3.89|), which can be studied in a way similar 



to Eq.(3.71). By repeating a discussion completely analogous to the one 
foll owing Eq.(3.71), we find that \e,h{v) — e ft°©l admits the same bound 

so that also lim^oo £ 2 {eh(v) — e£°(6)) =0. In 



Eq.(3.90) as \ef(v) - e^(v 



conclusion, the finite volume corrections Eq.(3.86) are exactly zero. Putting 



all the contributions together, we find that the central charge c defined by 



Eq.(3.17) is independent of A and equal to 1/2, as desired. 
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4 Proof of the partition function inequality. 



In this section we prove Lemma [TJ We start by proving it in the case A = 
with generic couplings {J b }- In this case, using Eqs.(2.4)-(2.7), we immedi- 
ately find 

Z%_ ({ J b }; A £>L ) + Z\({ J b }- A e , L ) = 2Z° e _ e ({ J b }; A e>L ) + 2Z o _ o ({J b }; A e , L ) 

(4.1) 

and 

Zl_{{J b }-K^ L ) + Z° + ({J b };A etL ) + Z°_({J b } ] A e<L ) = 3Z° e _ e ({ J b }; A l>L ) + 
+Z e _ o ({ J b }; A e>L ) + Z o _ e ({J b }; At >L ) + Z o _ o ({J b }; A ijL ) . (4.2) 

Remember that the four partition functions Z®_ e ({J b }; A^l), Z®_ ({ J b }; A^l), 
Zg_ e ({ J b }; A^l) and Z®_ ({ J b }; A^l) are all positive, see the definitions Eq. 



(2.8) and following lines. Therefore, the right hand side of Eq.(4.1) is > 



which proves Eq.(3.12) for A = and bond-dependent couplings. Moreover, 



the r.h.s. of Eq.(4.2) is bounded from above by 



3Z° e _ e ({ J b }; A i>L ) + 3Z e _ o ({ J b }; A lL ) + 3Z o _ e ({ J b }; A e , L ) + 3Z«_ {{ J b }; A £ , L ) 
= 3Z°({J b y,A. L ) , (4.3) 



where we used Eq.(2.9). This proves the lower bound in Eq.(3.11) for A = 



and bond-dependent couplings. The upper bound in Eq.(3.11) for this case 



is proved analogously: it is enough to observe that the r.h.s. of Eq.(4.2) can 
be bounded from below by 

Z° e _ e ({ J b }; A e>L ) + Zl ({J b }; A £>L ) + Z° _ e ({J b }; A e>L ) + Z o _ o ({J b }; A,, x ) 
= Z°({J b };A e>L ), (4.4) 

which leads to the desired bound. 

Let us now turn to the interacting case. The key issue is to obtain a rep- 
resentation of the four Grassmann partition functions Z ++ (Ai )L ), Z + _(A^ L ), 

-2L + (A^l), Z {A^l) in terms of (positive!) multipolygons partition sums. 

Such a representation is implicitly derived in [HJ Section II. A]. Let us make 
it explicit here. The starting point is the representation for the partition 
function in terms of disconnected polymers (see [TU Eqs.(2.8)-(2.9)-(2.10)]): 



ncosh 2 (^(x-y))j^e^ 



{*,y} 



E E ^({7i,---,7n})n^(7) 

™>0 {7i,...,7n}CA fji spgf 



(4.5) 
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where 

*(t)= e [n tanh (^ A ^)][ n *>] ( 4 - 6 ) 



and: 



5 connected: SGS bgbl(5) 
7(5)=7 



6 indicates a nearest neighbor bond and a^, is the bond spin, i.e. the 
product of the two spins at the vertices of b; 

S = {Si, . . . , S m } is a set of strings, where each string Si is the union 
of the bonds in a finite lattice path that can be either horizontal, or 
vertical, or "corner-like", as in Figure [l] moreover, v$ '■— f(x — y), 
where x and y are the first and last points connected by the path S on 
Ai t L. We say that a set of strings S = {Si, . . . , S m } is connected if, given 
1 < io < jo < m , we can find a sequence (Si , S^, . . . , Si p = Sj ) such 
that 5*4, PI Si l+1 7^ 0. From a graphical point of view, every connected 
component S corresponds in a non-unique way to a polymer 7(5), i.e., 
a connected set of bonds. It is helpful to color the bonds in 7(<S) black 
or gray, depending on whether the given bond belongs to an odd or 
even number of strings in S, and denote the set of bonds thus colored 
black by bl (5). See [H Fig.4]. 

the sum 7n } * n ^ ne r.h.s. runs over sets of polymers {71, . . . , 

such that each polymer is contained in A^ ^. Moreover, the function 
(^({Tl, . . . ,7„}) implements the hard core condition, that is if is equal 
to 1 if none of the polymers overlap, and otherwise (here two polymers 
overlap if and only if they have at least one bond in common); the term 
with n = should be interpreted as 1. 



Now, in order to obtain the Grassmann representation Eq.(3.1) one can 



proceed as explained in [T4"l Section II. A]: the idea is simply to rewrite 



Eq.(4.5) as 



= [ II cosh 2 A^(x -y))]-£ J2 ■ ■ ■ > In}) ■ 

{x,y} n>0 {7i,...,7„}CA4i 

n 

e [nn tanh (^)]E n ^ 



Si,...,S„ connected: i=l S£Si q_ bebl(>S) 

7(«5<)=7i 

where S := (S u ...,S n ), bl(<S) = Uf =l U(Si) and 

e n ^^=[ n -Ay^y-^y^ ^ 

a; beb\(S) 6ebl(5) P 
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Figure 1: The four possible types of "corner like" strings 



and then to re-express Z°({ Jb}', A^l) as a sum of Grassmann integrals, via 
Eqs. plQ] ) and Q. This leads to Eq.Q, with 



Z a (A e , L ) 

■ £ 



iSi,...,iS n connected: i=l S£c5; 
7(5i)=7i 



[ n cos^ (!/tt«(x - y))] • x; e ^> ■ ■ ■ > ^» 

{x,y} n>0 {7i,...,7„}CA £iL 

n 1 

II If t«nh(l/3X» s )] [ J] 

bebl(S) ^ 



2£({J 6 };A^) 



J h =.l 



(4.9) 



which is equivalent to jlU Eq.(2.11)] that, if further manipulated, implies 



the representation Eq.(3.2). On the other hand, if we plug Eqs.(2.4)-(2.7) 



into Eq.( 4.9[ ), we immediately get 



Z++(A^ L ) = Z e ^ e (Ae :L ) — Z e ^ (Ae :L ) — Z ^ e {A^ L ) — Z - (A^ L ) , 
<Z+_(A /jL ) = Z e - e (A tjL ) + Z e _ (A £iL ) - Z _ e (A £iL ) + Z _ (A 4L ) , 

Z-+(A EiL ) = Z e _ e (A /)L ) - 2 e -o(A^ L ) + £ - e (A £jL ) + Z _ (A /)L ) , 

Z__(A^ L ) = Z e ^ e (Ai >L ) + Z e _ (A £jL ) + Z _ e (A £iL ) - Z _ (A| )X ) , 

(4.10) 
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with 



Z e ^ e (A l>L ) = [ J] cosh 2 Q/3A*;(x-y))] • £ £ ^({ft, . . . , %}) 



{x,y} 



™>0 {7i,...,7»}CA^z, 

1 a 1 



e [n n tanhd^)] [ n ^k.cwjm^ ^.n 



iSi,...,iS n connected: i=l SGiS,: 
7(<S;)=7i 



bSbl(S) 



and analogously for the three other partition functions, with the label e — e 
replaced at both sides by e — o, o — e, — 0, respectively. If we use the 



definition Eq.rt2.8b, we can rewrite the last factor in Eq.(4.11) as: 



n \U Z ' UW ^ ] 

b&l{S) r 



\J h =J 



2 ^|bi(5), [ -Q cosh(/3 Jb)] ^ Yl Yl m , (4.12) 



where is equal either to t, if b E~ Bl(iS) = Us. e sbl(«St), or to 1/t, if 

b E Bl(5). An equivalent way of rewriting Eq.( |4.12 ) is 



n 

6ebl(5) 



ld_ 
PdJ b 



Z e °_ e ({J 6 };A,, L ) 



j,=j 



t Mm\ z l e ({J b (S)y,A i>L ), (4.13) 



where J&(£) is equal either to J, if b Bl(5) = Us ie sbl(«Si), or to (3 1 arctanh(l/t), 
if b E Bl(5). Plugging this back into Eq.( 4.1lj ) gives 



Z e _ e (A e , L ) = [ H cosh 2 Q/3A*;(x-y))] £ <p({% ■ . . , In}) 



{x,y} 

n 



e [nn tanii (w 



™>0 {7i,...,7„}CA^£ 

£ |bi(5)l^o_ e({ J 6( 5 )};A4i) 



(4.14) 



S\,...,S n connected: i=l S(ESi 
7(<%)=7i 



and analogous formulas are valid for the three other partition functions, with 
the label e — e replaced at both sides by e — o, o — e, o — o, respectively. 
Now, the key observation is that if Xvs > for all S, then all the factors 



appearing in Eq.(4.14) (as well as in its analogues with e — e replaced by 



e — o, o — e, — 0) are non-negative. Therefore, if we insert Eq.(4.14) and its 



analogues with the label e — e replaced by e — o, o — e, o — o, into Eq.(4.10) 



and if we use the known bounds on the partition functions at A = with 
bond-dependent couplings, proved at the beginning of this section, namely 

Z\{ J b (S)}; A e>L ) < Z°__({J b (S)}; A^ L )+Z\{{ J b (S)}; A^ L )+Z° + _({ J 6 («S)}; A^ L ) 
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Z°_({J 6 (5)};A^)+Z° + ({J b (5)};A <iL )+Z°_({J 6 (5)};A,, L ) < 3Z°{{J b {S)}; A e>L ) 
Z°_ + ({J b (S)};A eiL ) + Z° + _({J b (S)};A e>L ) > 



then we finally obtain the desired estimates Eqs.(3.11 )-(3.12). This concludes 



the proof of Lemma [T] and, therefore, of Theorem LI 



A The contribution to the pressure from the 
ratio of the Grassmann partition functions 

A.l The non- interacting case 

In this section we prove the vanishing of the limit in the second line of 



Eq.(2.17) 



lim lim — log 



1 Z°_ + (A e>L ) Z° + _(A ijL ) 
2 [ + Z°_(A £>L ) + 



/3=/3 c J 







(A.l) 



To this purpose, we rewrite the partition functions in Eq.(2.15) as 

i-i l-i 



2) TI 11(2 



M2r+1 — Z2n+l 



Z 2n+l) 



1 1/2 



r=0 n=0 

e-i l-i 



(V2Y L H \H(2a 2r+1 -z 2n -z£) 



r=0 n=0 

e-i l-i 



1/2 



(A.2) 



(V2Y L H l[(2a 2 . 



2r z 2n+l Z 2n+1/ 



1/2 



r=0 n=0 



where a p := 2 — cos(|p) and z p := e l i p . The expressions in square brackets 
in the r.h.s of these equations can be further rewritten and put in the form 
used by [10J. Consider first the expression in square brackets appearing in 

the definition of Z°__{Ai^L) . Note that z 2n +i are the L roots of —1 and, 



therefore, n^ol- 2 ~ z 2n+i) = z L + I, \fz E C. In particular, n«=o z 2n+i 



L-l 



[-1) L , so that 

L-l 



L-l 



J^(2a 2r +1 — ^2n+l — Z 2n+l) = Jlfen+l ~ ^ a 2r+l z 2n+l + 1) = 
n=0 n=0 
L-l 

= Y[( Z 2n+l ~ a tr+l)( Z 2n+l ~ a 2r+ l) = [(4r+l) L + X ] ' [i a t+l) L + X ] 
n=0 
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where a p = a p ± */a 2 — 1. Choosing 7 P > in such a way that cosli7 P = a p , 
we can further simplify this into 

L-l 



n=0 



-1 



«2r+l — 2 2n+l ~ Z 2n+1 



e L 72r+1 + !)( e -^+i + !) = 4 cos h 2 



(A.3) 



Plugging this back into the definition of Z°__(A^i) 



e-i 



P=Pc 



IL 



\\ 2 cosh 



P=Pc 



gives 



r=0 



which is the same as [TUl Eq.(2.3)]. Actually, by exchanging the roles of 
and L we can rewrite Z°__(Ai^)\ in two equivalent ways: 



P=p c 



P=Pc 



L-l 



2r+l 



CL 



Y[ 2 cosh 



/2n+l 



r=0 



n=0 



(A.4) 

where 7 P > is defined by the condition that cosh7 P = 2 — cos(-^p). Pro- 
ceeding exactly in the same way for Z°_, (A e A , Z\_ (A^ A gives: 

' P=Pc ' \P=Pc 



P=Pc 



(V2Y L Y[2smh 



Lj, 



L-l 



2r+l 



l2v 



r=0 

e-i 



zi_(A ijL ) _ =( V ^r i n 2cosh ^? : =(^) £L n 2sinh 



n=0 
L-l 



2 

2n+l 



P=Pc 



r=0 



n=0 



By using these formulas we rewrite the l.h.s. of Eq.(A.l) as 

L-l 



lim lim — log 

I— >oo L— >oo .L 



1 

L2 



(1+ntmh ^ + ntm h^±i) 



r=0 



ra=0 



(A.5) 



Since < tanhx < 1 for x > 0, it is apparent that the argument of the 
logarithm in this equation is positive and smaller than 3/2, which implies 



that the limit in Eq.(A.5) is zero. 



A. 2 The interacting case 



In order to bound the ratios 



\Z-+(A tiL )\ 
Z—(A e , L ) 



and ^ + ^ e ' L l^ appearing in Eq. (13.151) 



we compute Z_ + (Ae yL ) and Z + _(A e>L ) by a renormalization group construc- 
tion analogous to the one used to analyze Z (A^l). Everything is the same, 
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with a few obvious changes induced by the different boundary conditions. 
The important fact that makes the renormalization group construction of 
Z_ + (Ai tL ) and Z + _(A^ L ) possible, is that h* is finite for both: in the first case 

h* = \}og 2 ('x / £) \ , as for Z (A^l), while in the second h* = |_log 2 (7r/L) J . 

Of course, the bulk contributions to the free energy are the same for all these 
partition functions, so that 



\Z-+(Ai, L )\ < \Z1 + {All 



\Z + -{^,l)\ 



< 



2L_(A<, L ) " YZ*__{A LL 



Z°-(All) 
Z + -(Ae, L ) 



JLR 2 (X) 



(A.6) 
(A.7) 



where Z°^_(A 



LL) 



Z°_ + (A e 



Z ( 1_(A£ i l)\ /3=i3 are the non interact- 



ed 



ing partition functions at criticality defined in Eqs.(2.14)-(2.15), and Ri, R 2 , 
according to be renormalization group analysis, are bounded by c|A|£~ 2 ~ 9 
for some c > and < 6' < 1. As proved in Appendix |Al the ratios 



Z°+(^.l) 

therefore, 



and 



.(A,i) 



are positive and smaller than 1 and, 



\Z-+(At 
Z_4A 



+ 



LL, 



Z_4A 



< e 



C\X\i~ 



LL) 



(A.8) 



for a suitable constant C > 0. Plugging this back into the l.h.s. of Eq.(3.15) 
gives 



£ / 
lim lim — log ( 1 + 



|Z_ + (A,, L )| | |Z + _(A <|L 

Z-{^L) 



, < lim lim — ClAl - , a . , 

Z__(A^ L ) J ~ l^ooL^ooL ' 

(A.9) 



which proves Eq.(3.15). 
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